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Abstract 

A  two-dimensional  transient  simulation  model  for  a  cascaded  (10  cascades)  micro-tubular  solid  oxide  fuel  cell  (SOFC)  and  a  common  micro¬ 
tubular  SOFC  with  a  preheater  tub  was  developed  to  calculate  the  distribution  of  gas  species,  the  local  current,  voltage  and  temperature.  The 
goal  was  to  compare  these  two  different  SOFC  designs  under  steady  state  and  load  change  conditions  with  respect  to  power  density.  Steady  state 
simulation  results  have  shown  that  the  cascaded  cell  concept  is  able  to  operate  at  higher  average  cell  voltages  in  principle.  Compared  to  the  cascaded 
cell  concept  the  average  ohmic  anodic  and  cathodic  resistance  of  the  common  cell  has  to  be  84  times  lower  to  reach  the  same  average  power  density 
at  an  average  voltage  of  0.7  V.  The  cascaded  cell  concept  gives  the  opportunity  to  operate  at  a  higher  average  voltage  than  the  lowest  Nernst  voltage 
offers.  This  could  be  interesting  to  realize  high  fuel  utilization  (>85%).  Transient  simulation  results  have  shown  that  both  cell  concepts  respond 
within  10  s  on  an  electrical  load  change.  The  cascaded  cell  concept  has  shown  a  lower  temporary  voltage  drop  than  the  common  cell  concept.  This 
paper  describes  the  theory,  solution  techniques  and  results. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  micro-tubular  SOFC  has  advantages  due  to  rapid  start-up 

[1] ,  load  change  behavior  and  the  volume  related  power  density 

[2]  because  the  cell  area  per  volume  increases  by  decreasing 
the  cell  diameter.  The  thermal  integration  of  heat  sinks  is  an 
important  issue  to  simplify  SOFC  systems.  The  large  scaled 
tubular  SOFC  developed  by  Westinghouse  consists  of  an  air 
preheater  tube  to  preheat  the  air  directly  by  the  generated  heat. 
This  concept  gives  the  further  opportunity  to  implement  heat 
sinks  like  a  stack-integrated  prereformer.  Thereby,  the  excess 
air  demand  is  reduced.  Thermodynamic  and  economic  analyses 
have  shown  that  the  low  excess  air  demand  and  stack-integrated 
heat  sinks  are  a  prerequisite  to  design  low  costs  and  efficiently 
SOFC  systems  [3].  Therefore,  the  air  preheater  tube  is  adopted 
for  this  investigation  of  micro-tubular  SOFCs  (Fig.  1). 

The  tubular  SOFC  of  Siemens- Westinghouse  has  compared 
to  the  planar  SOFC  the  disadvantage  that  the  ohmic  losses  are 
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much  higher  due  to  the  longer  ways  of  electrical  current.  The 
micro-tubular  concept  developed  by  Adelan  Ltd.,  UK,  uses  wires 
to  collect  the  current.  This  reduces  the  internal  ways  of  electrical 
current  in  comparison  to  the  tubular  concept  of  Westinghouse. 
The  planar  concept  provides  the  largest  area  for  the  current  col¬ 
lector  but  it  shows  disadvantages  in  start-up  time  and  mechanical 
stability.  Flowever,  the  ohmic  resistance  of  a  current  collector  is 
an  important  issue,  which  could  have  a  strong  influence  on  power 
density. 

The  local  cell  voltage  is,  from  a  thermodynamic  point  of 
view,  limited  by  the  local  Nernst  voltage.  This  means  that  the 
maximum  voltage  of  common  single  cells  with  equipotential 
electrodes  could  reach  the  lowest  Nernst  voltage  (highest  fuel 
utilization).  Cascaded  concepts  (Fig.  2)  give  the  opportunity  of 
realizing  different  voltages  at  the  cascades.  These  voltages  are 
added  up  to  the  all  over  cell  voltage.  If  the  cascaded  average  cell 
voltage  has  the  same  value  as  the  single  cell  voltage,  the  entropy 
production  must  also  have  the  same  value  as  well. 

Cascaded  concepts  have  the  advantage  of  reaching  higher 
average  voltages  than  the  lowest  Nernst  voltage  in  the  gas  bulk. 
On  the  other  hand  current  efforts  are  focussing  on  average  volt¬ 
ages  which  where  lower  than  the  lowest  Nernst  voltages  at  a  fuel 
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ch4,  h2o,  H2,  CO,  co2 


Fig.  L  Micro-tubular  SOFC. 


e' 


CH4,  H20,  H2,  CO,  C02 
Fig.  2.  Cascaded  micro-tubular  SOFC. 


utilization  of  85%.  This  leads  to  higher  average  power  densities. 
In  case  of  an  increasing  fuel  utilization  the  cascaded  concepts 
could  be  an  interesting  solution  to  reach  layout  voltages  which 
where  higher  than  the  lowest  Nernst  voltage.  Compared  to  com¬ 
mon  single  cells,  cascaded  concepts  are  reducing  the  all  over 
electrical  current  by  the  number  of  cascades.  This  could  lead 
to  lower  ohmic  polarization  losses  than  the  common  single  cell 
concept. 


Fig.  3.  Schematics  of  heat  transfer. 


constant  and  equal  in  each  direction. 

1  dT  d2T  d2T 


div  grad  T  — 


r  dr  dr 2  dz2 


(2) 


The  differential  equation  is  applied  in  axial  and  radial  direc¬ 
tion  of  the  cell.  The  temperature  profile  in  the  circumferential 
direction  is  assumed  to  be  homogenous. 

There  are  heat  sources  due  to  the  reaction  enthalpy  of  the 
reforming,  the  shift  reaction  (gref+shift)  and  the  reversible  heat 
(<?ox,rev)  of  fuel  oxidation,  respectively. 

1 

<7reac  =  “ - .  .  (A//ref  +  A  shift  +  A  Q  ox,  rev  )  (3) 

A  r  A  A 

This  heat  is  equal  to  the  difference  of  the  all  over  reaction 
enthalpy  and  the  reversible  work  of  fuel  oxidation. 

1 

<7 reac  —  ,  “  “( A//reac(7"=const)  fF0x,rev)  (4) 

A  r  A  A 


2.  Model 

2.1.  Energy  balance 

For  the  determination  of  the  cell’s  electrical  performance  it 
is  necessary  to  know  the  cell’s  temperature  distribution  because 
the  loss  mechanisms  are  strongly  temperature  sensitive.  Hall 
and  Hall  and  Colclaser  [4]  proposed  a  one-dimensional  transient 
model  of  a  tubular  SOFC,  which  uses  Eq.  (1)  to  calculate  the 
temperature  distribution  of  the  solid  ceramic  of  the  SOFC.  Heat 
transfer  (Fig.  3)  occurs  through  thermal  conductivity,  convection 
( qa )  and  radiation  (qrM\).  The  conductive  heat  transfer  in  the  solid 
material  of  the  cell  can  be  calculated  by  solving  the  local  energy 
balance. 

dT  q 

—  —  a  div  grad  T  H -  (1) 

dr  pep 

Here  T  is  the  temperature,  r  the  time,  p  the  solid  mass  density, 
cp  the  specific  heat,  a  the  thermal  conductivity  coefficient  of  the 
solid  ceramic  of  the  SOFC  and  q  is  the  volume  specific  heat 
source.  It  is  assumed  that  the  thermal  conductivity  coefficient  is 


The  oxidation  of  hydrogen  and  carbon  monoxide  are  assumed 
to  be  the  only  electrochemical  reactions.  It  can  be  shown  that 
the  reversible  work  of  these  two  reactions  have  the  same  value 
if  the  shift  reaction  is  fulfilling  the  equilibrium  condition.  Inves¬ 
tigations  of  the  reforming  reaction  of  methane  at  Ni-cermet 
anodes  [10]  have  shown  that  the  shift  reaction  is  near  to  equilib¬ 
rium  condition.  The  different  mechanism  losses  of  the  hydrogen 
and  carbon  monoxide  oxidation  were  neglected  to  determine 
the  reversible  work  from  only  one  of  these  reactions.  The 
following  Eq.  (5)  uses  the  reversible  work  of  the  hydrogen 
oxidation 


1 

<7reac  —  “  ~ 7 

Ar  AA 


(  Tsr  \ 

5>out 

A  bH?+  cpJ(T)dT 

i 

\  T0  ) 

(  Tsr  ) 

A bH°+  /  cp  i(T)dT 

—  ^H2,ox^ox,rev 

V  i  J 

(5) 
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which  can  be  determined  by  the  free  reaction  enthalpy  of  the 
hydrogen  oxidation. 


Wox,rev  =[(  A  rG{T)  +  TR  In 

J  \  pn2(uf)^/po2(uf) 

uf,  1  N  ' 


Convection  occurs  between  solid  surfaces  and  the  gas 
streams. 


Qa  — 


a(Nu) 
A  r 


(//solid  fluid)  — 


1 

Ar  AA 


A //fluid 


(7) 


where  A  r  the  radial  segment  and  AA  the  area  segment  perpen¬ 
dicular  to  flow  direction  of  fuel  and  air. 

It  is  assumed  that  radiation  transfers  only  occur  between  the 
cathode  and  air  preheater  tube. 


f/  rad  — 


C  Ca,ft 
Ar 


(8) 


where  x  is  the  molar  fraction  of  reacting  species,  A  rH  is  the 
reaction  enthalpy  and  A  ,  S  is  the  reaction  entropy.  The  inves¬ 
tigations  of  the  shift  reaction  at  Ni-cermet  anodes  have  shown 
[  10,1 1]  that  the  conversion  rates  are  always  near  to  equilibrium 
condition. 

The  reforming  reaction  is  a  kinetic  controlled  reaction.  Some 
investigations  [9]  at  Ni-cermet  anodes  shown  that  boundary  dif¬ 
fusion  and  pore  diffusion  can  be  neglected  below  700  °C  and 
boundary  diffusion  can  be  neglected  between  800  and  900  °C 
[  12].  In  this  study,  both  transport  processes  were  neglected.  The 
conversion  of  methane  is  given  [10]  as 

«CH4,0  -  «CH4  =  A  AKat  hf{  PCYU,  Pu20  (1 1) 

where  AAicat  is  the  catalytic  area  calculated  by  the  geometry  of 
the  anode,  kt g  the  effective  rate  coefficient  and  p  the  partial  pres¬ 
sure  of  gas  species.  The  conversion  rate  coefficient  is  determined 
by  the  Arrhenius  equation 

keff  =  k0  exp  (!2) 


Here  is  C  the  radiation  exchange  factor  between  cathode  and  the 
preheater  tube. 

The  heat  sources  due  to  the  polarization  effects  (qn)  are 
expressed  by 

fji  A  A,  R(i  AA,)2  8R$i 2  A  A;  ^ 

^  ~  Ar  AA  ~~  Ar  AA  ~  Ar  AA  (9) 

where  r)  is  the  overpotential,  i  the  current  density,  AA,  the  area 
segment  in  the  equivalent  electrical  current  flow  direction,  R 
the  equivalent  ohmic  resistance  in  £2,  R$  the  equivalent  specific 
ohmic  resistance  in  £2  m  and  8  is  the  thickness  of  equivalent 
ohmic  resistance. 

2.2.  Mass  balance 


where  ko  is  the  pre-exponential  coefficient  and  E,\a  the  activa¬ 
tion  energy.  The  boundary  of  conversion  rate  is  given  by  the 
thermodynamic  equilibrium: 


*h2*co  (P°Y  (  ArH  —  T  ArS 

-*ch4  V'II20  V  P  J  exp  V  RT 


(13) 


The  local  conversion  of  the  hydrogen  oxidation  is  determined 
by  Faradays  law  as  follows: 


«h2,0  —  «h2 


«O2,0  -  « o2 


i  AA 
2  F 
i  AA 
4  F 


(14) 

(15) 


For  the  calculation  the  fluid  channel  is  divided  into  slices 
in  the  same  distances  as  assumed  for  the  cell.  For  the  sim¬ 
plification  the  gas  composition  is  assumed  to  be  homogenous 
distributed  in  each  gas  slice.  The  chemical  species  CH4,  Ho, 
H2O,  CO,  CO2  for  fuel  and  O2,  N2  for  oxidant  are  consid¬ 
ered  in  this  study.  It  is  assumed  that  all  the  gases  are  ideal 
and  incompressible  fluids.  The  considered  chemical  reactions 
are: 

CH4  +  H2O  -o-  3H2  +  CO  Steam  reforming 
CO  +  H2O  -o-  CO2  +  H2  Shift  reaction 
H2  +  O.5O2**  H2O  Hydrogen  oxidation 


The  all  over  fuel  utilization,  uf,  can  be  calculated  by  the  cell 
current  /,  number  of  cascades  ncas  and  inlet  gas  composition 

I  Ur4s 

uf= - — -  (16) 

C(8«ch4,o  +  2«co,o  +  2«h2.o) 

In  case  of  steady  state  calculations  the  inlet  mass  flow  is  deter¬ 
mined  iteratively  to  reach  a  requested  fuel  utilization. 

2.3.  Electrochemical  approaches 

Polarization  is  a  common  parameter  in  the  analysis  of  fuel 
cell  performance.  To  calculate  the  relationship  between  local 
voltage  I/celi.z.  and  local  current  density  the  following  general 
expression  can  be  used: 


In  this  study,  the  oxidation  of  carbon  monoxide  is  neglected.  The 
reversible  work  of  the  hydrogen  oxidation  and  carbon  monoxide 
oxidation  have  the  same  value  if  the  shift  reaction  fulfills  the 
equilibrium  condition 


v  1 1 2  xCq2  _  (  ArH-TArS 

xco  xh2o  V  RT 


(10) 


Ucell.L  —  C  NJcrnst  f/Act.  An  f/Act,Ca 

-  ^Diff.An  -  ^Diff.Ca  —  r]S2. Elt  (17) 

The  overpotential  can  be  divided  in  activation,  diffusion 
(concentration)  and  ohmic  overpotential.  The  Nernst  voltage  is 
considered  to  be  the  driving  force,  which  is  a  function  of  the 
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reversible  work  of  the  hydrogen  oxidation  (Eq.  (6)). 


jj  _  Wox.rev 

t^Nemst  =  7  _ 

Auf2F 


(18) 


For  calculating  the  resulting  voltage  of  one  cascade  further 
ohmic  resistances  (Rn, An.  Rn, Ca  as  shown  in  Fig.  4)  have  to  be 
taken  into  account.  The  voltage  drop  over  these  ohmic  resis¬ 
tances  is  determined  iteratively  by  using  Kirchhoff’s  law  and 
the  local  voltages  t/ceii,L- 


where  FAct.An  is  the  activation  energy,  yAn  an  empirical  pre¬ 
exponential  factor,  p\\2  the  partial  pressure  of  hydrogen,  pn20 
the  partial  pressure  of  water  and  po.An  is  the  reference  pressure 
at  the  anode. 

The  cathodic  exchange  current  density  of  FSM  cathodes  is 
expressed  as  [13]: 


'0,Ca  =  KCa 


PO 2_\ 
PO.  Ka  / 


0.25 


exp 


f-Act.Ca  \ 

RT  ) 


(22) 


2.4.1.  Concentration  overpotential 

The  concentration  overpotential  is  a  Active  electrical  potential 
drop  caused  by  the  concentration  difference  in  the  gas  phase. 
Diffusion  through  the  porous  material  is  typically  described  by 
binary  molecular  diffusion  and  Knudsen  diffusion.  In  this  study 
both  types  of  diffusion  are  taken  into  account  to  estimate  the 
concentration  overpotential.  It  is  assumed  that  the  concentration 
gradient  along  the  pore  is  much  higher  than  the  concentration 
gradient  inside  the  surface  boundary.  Only  the  concentration 
gradient  in  the  porous  electrodes  is  taken  into  account. 

For  the  ideal  gas  model  the  anodic  concentration  overpoten¬ 
tial,  which  is  proportional  to  the  irreversible  entropy  production 
can  be  written  as 

_  RT  ,  ( Fh2  Ph2o\  _  T  a  c 
“  ZF  ln  ( ,4  p'H20 )~zF  AS”-A"  <23) 


2.4.  Activation  overpotential 


The  activation  of  a  chemical  reaction  involves  energy  barriers, 
which  must  be  overcome  by  the  reacting  species.  These  energy 
barriers  are  called  the  activation  energy  and  result  in  activation 
or  charge  transfer  polarization.  The  activation  polarization  of  an 
electrochemical  reaction  is  usually  expressed  by  the  well-known 
Butler- Volmer  equation 


i  =  i  o 


exp 


uzFriAcA 
RT  ) 


(1  -  &)zFf)Act  V 
RT  )_ 


(19) 


and  for  the  cathodic  concentration  overpotential  as  follows: 


RT 

hDiff.Ca  =  t—  In 
2  zF 


—  -  ASin-  Ca 
2  zF 


(24) 


where  p°  is  the  partial  pressure  in  the  gas  bulk,  p *  the  partial  pres¬ 
sure  at  the  end  of  the  pores.  The  concentration  profiles  inside  the 
pores  can  be  estimated  by  using  the  Fick’s  law  (Eqs.  (25)— (27)). 
The  molar  fraction  x  can  be  expressed  as 


*  _  o  _  RTi 8 An 

X"2  ~  Xli 2  p2FDetfiAn 


(25) 


where  a  is  the  transfer  coefficient  and  z'q  the  exchange  current 
density.  The  transfer  coefficient,  which  is  approximately  0.5 
for  fuel  cells,  is  considered  to  be  the  fraction  of  the  change 
in  polarization  that  leads  to  a  change  in  the  reaction  rate  con¬ 
stant  and  its  value.  The  exchange  current  density  is  the  cathodic 
and  anodic  electrode  reaction  rate  at  the  equilibrium  potential.  A 
high  exchange  current  density  means  that  good  fuel  cell  perfor¬ 
mance  can  be  expected.  With  a  =  0.5  (Eq.  (19))  the  activation 
polarization  can  be  calculated  as  follows: 

RT  (  i  \ 

=  0A?  “p  UJ  <20) 

For  the  anodic  exchange  current  density  of  the  hydrogen  oxi¬ 
dation  at  Ni-YSZ  anodes,  Mogensen  and  Lindegaard  [5]  has 
proposed  Eq.  (21). 


(-™D)  ('£HSo')exp('_£AelA,'l 
\P0.AnJ  \P0,AnJ  \  RT  J 


(21) 


'  h2o 


'h2o 


RTi  8 An 
p2FDeff  An 


for  the  equimolar  diffusion  of  hydrogen  and  water  and 


Xq2  =  1  +  (*o2  -  1)  exp 


/  RTiScn  \ 

Vp4FDeff,Ca/ 


(26) 


(27) 


for  the  self-diffusion  of  oxygen  in  nitrogen.  ^An  is  the  thickness 
of  the  anode,  8q a  the  thickness  of  the  cathode  and  /.)err  the  effec¬ 
tive  diffusion  coefficient.  In  this  model,  the  effective  diffusion 
coefficient  takes  the  Knudsen  diffusion  as  well  as  the  ordinary 
diffusion  at  the  cathode  and  the  binary  diffusion  at  the  anode 
into  account. 

The  binary  diffusion  coefficient  between  gas  species  A  and  B 
in  free  space  is  given  by  the  Chapman-Enskog  theory 

/j  1  \’/2  r3/2 

Dab  =  1.8583  x  1(T6  —  +  — 

\Ma  Mb  J 


*0,An  =  yAn 


po\BT2D 


(28) 
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where  a  is  the  collision  diameter  (Lennard-Jones  length)  and 
Q[)  is  the  collision  integral  based  on  the  Lennard-Jones  poten¬ 
tial.  Neufield  et  al.  [6]  has  proposed  the  following  approach 

^  1.06036  0.19300 

°  ~  (kT/sAB)015610  +  exp(0.47635(A-77eAB)) 

1.03587  ;  1.76474 

+  exp(1.52996(ytr/eAB))  +  exp{3.894U(kT/sAB)) 

(29) 


to  calculate  the  collision  integral  by  the  Lennard-Jones  energy  e, 
which  is  a  function  of  the  critical  gas  properties  and  Boltzmann 
constant  k. 

The  Knudsen  diffusion  coefficient  for  a  species  A  is  defined 
as 


2J?pori>A 
i-'Kn,  A  =  - “ - 


(30) 


where  Rpor  is  the  pore  radius  and  v  is  the  average  velocity  of  the 
gas  molecule.  Both  ordinary  diffusion  and  Knudsen  diffusion 
can  be  combined  within  the  effective  diffusion  coefficient  which 
can  be  written  as: 


AjIT.Por.A  0  \DaB  Dan,  A 

The  tortuosity  rp  and  porosity  9  are  used  to  account  the  geom¬ 
etry  of  the  pores  and  relation  between  solid  space  and  pore  space 
of  the  electrode. 

If  we  neglect  the  presence  of  other  species  except  hydrogen 
and  water  inside  the  anodic  pores,  the  anodic  effective  diffusion 
coefficient  can  be  expressed  for  the  ideal  gas  mixture  as  [7]: 


Arff.An  =  (  1  H2C>N)  T>eff,Por,H2  +  ( Deff,Por,H20  (32) 
V  Fges  J  V  Pges  / 

Inside  the  cathodic  pores  the  cathodic  effective  diffusion 
coefficient  can  be  formulated  as 


2.5.  Solution  technique 

The  model  includes  the  calculation  of  chemical  reactions 
like  the  steam  reforming,  shift  reaction,  electrochemical  reac¬ 
tions  and  thermal  heat  transfer.  The  thermal  part  of  the  model 
accepts  the  cell  electrical  output  and  the  cell  heat  sources  from 
the  electrochemical  part  and  calculates  new  cell  temperature. 
Overpotential  and  Kirchhoff’s  law  are  used  to  determine  the 
local  current  density  and  local  cell  voltage.  The  local  cell  cur¬ 
rents  are  summed  up  to  the  total  cell  current.  The  cell  voltage  is 
determined  by  the  sum  of  cascade  voltages.  For  the  calculation 
of  steady  state  conditions  the  calculated  cell  voltage  is  compared 
with  the  requested  cell  voltage.  Cell  heat  sources  are  calculated 
from  reversible  entropy  change  due  to  the  oxidation  of  hydro¬ 
gen,  the  local  overpotential,  the  heat  transfer  by  radiation  and  the 
heat  transfer  by  convection,  which  includes  the  sensible  enthalpy 
changes  in  the  gaseous  phases.  These  heat  sources  are  applied  to 
the  solid  nodes  of  the  thermal  SOFC  model  to  calculate  the  new 
cell  temperature  by  solving  Eq.  (1).  The  transient  model  uses 
the  Cranck-Nicolson  discretization  (axial:  51  nodes;  radial:  7 
nodes),  which  has  the  advantage  that  no  stability  boundaries 
have  to  be  considered.  The  generated  2D-tridiagonal  matrix  is 
solved  iteratively  by  the  Gauss-Seidel  method.  Fig.  5  shows  a 
flow  chart  of  this  process.  This  model  can  also  calculate  steady 
state  operation  points  by  running  transient  simulations  until  the 
transient  decay  to  a  level  below  a  specified  amount. 

The  model  is  used  to  determine  the  transient  response  as  a 
consequence  of  the  electrical  load  changes.  The  cell  is  subjected 
to  a  step  of  increasing  current  and  a  time  dependent  increment 
of  input  molar  gas  flows 

h{  t)  =  h(  r  =  0)  +  [h{x  =  oo)  —  h(r  —  0)](1  —  exp(— r/r  c)) 

(38) 

where  r  is  the  time  after  the  load  change  and  tc  is  the  time 
constant.  Hereby  the  input  gas  compositions  are  held  constant. 


Arff.Ca  =  7Jeff,Por,02 


(33)  3.  Results 


2.4.2.  Ohmic  overpotential 

Ohmic  losses  occur  according  to  the  flow  of  ions  in  the  elec¬ 
trolyte 

fin,  Eit  =  iRn,  Eit  (34) 

and  the  flow  of  electrons  through  the  electrode  materials. 

fin,  An  —  iRn,  An  (35) 

fin,  ca  =  iRn,  Ca  (36) 

The  ionic  conductivity  of  the  electrolyte  depends  strongly  on 
material,  temperature  T  and  thickness  8.  For  YSZ  electrolytes 
Achenbach  et  al.  [8]  has  proposed  the  following  dependency: 

1  £2  m 

Rq  Eit  =  -  (37) 

’  33,  400  exp(— 10,  300 K/  T) 

In  this  study,  the  electronic  resistances  of  the  anode  and  cath¬ 
ode  are  assumed  as  constant. 


The  current  and  temperature  distributions  at  the  steady  state 
and  transient  conditions  are  simulated  for  both  the  common  sin¬ 
gle  micro-tubular  cell  and  the  cascaded  micro-tubular  cell.  The 
initial  parameters  are  listed  in  Table  1 .  In  this  study,  temperatures 
of  the  fuel  and  the  oxidant  at  the  inlets  of  the  cell  are  assumed 
to  be  the  same. 

3.1.  Steady  state  characteristic 

The  temperature  distribution  of  the  common  micro-tubular 
cell  is  shown  in  Fig.  6.  In  comparison  to  the  results  of  large  scaled 
tubular  SOFC  reported  by  Siemens- Westinghouse  [  14]  the  qual¬ 
itatively  temperature  profile  of  the  micro-tubular  SOFC  shows 
similar  tendency.  However,  the  deviation  between  experimental 
data  and  simulation  data  indicates  that  the  future  comparison  to 
measured  data  of  micro-tubular  SOFC  is  necessary. 

The  molar  fraction  profile  along  the  flow  direction  (z- 
direction  from  fuel  inlet)  is  shown  in  Fig.  7. 
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Specify: 

2-D  geometry  (cartesian/cylindrical  coordinate  system) 
number  of  cascades 

discretization  method  (steady  state  /  transient) 
material  properties,  inlet  flow  conditions 
solid  temperaturprofile  for  transient  calculation  (x  =  Os) 
requested  voltage  for  steady  state  calculation 


I 


f 

Cranck-Nicolson 

\ 

V 

discretization 

/ 

X 


mass  balance 

-  reforming  reaction 

-  shift  reaction 

-  hydrogen  oxidation 


X 


electrochemical  balance 

local  Nemst  Volage 
-  local  current  density 
•  local  polarization 


energy  balance 

-  electrochemical  heat  sources 

-  heat  transfer  by  convection 

-  heat  transfer  by  radiation 

.  I 


Gauss-Seidel  method 
to  solve  the  2D-tridiagonaI 
matrix  for  the  temperature 
profile  in  the  solid  phase 


convergence 
yes 


outout: 

-  local  temperature 

-  local  voltage 

-  local  current 

\ 

-  gas  composition 

_ / 

Fig.  5.  Flowchart  of  solution  technique. 


In  the  first  section  of  the  fuel  cell,  the  hydrogen  fraction 
increases  due  to  the  endothermic  reforming  of  methane.  This 
has  the  effect  of  cell  cooling  at  the  fuel  inlet.  The  decreasing 
of  hydrogen  fraction  occurs  by  the  electrochemical  oxidation 
of  hydrogen,  which  is  proportional  to  the  local  current  density 
(Fig.  8).  Thereupon,  the  shift  reaction  produces  hydrogen  as  long 
as  carbon  monoxide  is  available. 

As  shown  in  Fig.  8,  the  local  Nernst  voltage  decreases  due 
to  this  drop  of  hydrogen  and  carbon  monoxide.  In  the  common 
micro-tubular  cell,  the  local  current  density  increases  along  the 
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Fig.  6.  Temperature  profile  of  the  common  micro-tubular  cell. 


Fig.  7.  Hydrogen  and  methane  fraction  of  the  common  cell. 


Fig.  8.  Local  voltage  and  current  density  of  the  common  cell. 
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Table  1 


Initial  parameters 


Number  of  cascades 

10 

Cell  length  [m] 

0.1 

Cell  inside  diameter  [m] 

2  x  10~3 

Thickness  anode  (Ni-YSZ)  [m] 

0.05  x  10~3 

Thickness  cathode  (Lao.6Sro.4Mn03)  [m] 

0.05  x  10~3 

Thickness  electrolyte  (YSZ)  [m] 

0.1  x  10“3 

Heat  capacity  ceramic  [Jkg-1  K-1] 

400 

Heat  conduction  coefficient  ceramic 

2 

[Wm_1  KA1] 

Density  ceramic  [kgm-3] 

6600 

Air  preheater  tube  diameter  [m] 

1  x  10~3 

Thickness  air  preheater  tube  [m] 

0.1  x  10-3 

Heat  conduction  coefficient  air  preheater 

6 

tube  [W  m_ 1  K-1] 

Pre-exponential  coefficient  for  exchange 

18.9  x  109 

current  density  anode  [Am  2] 

Pre-exponential  coefficient  for  exchange 

1034  x  109 

current  density  cathode  [Am  2] 

Activation  energy  anode  [kJ  mol-1] 

120 

Activation  energy  cathode  [kJ  mol-1] 

130 

Activation  energy  reference  [kJ  mol  - 1  ] 

150 

Ohmic  resistance  anode  cascaded  [£2  m] 

5.6  x  10~6 

Ohmic  restistance  cathode  cascaded  [£2  m] 

5.6  x  lO^6 

Radiation  exchange  factor  [W  m-2  K-4] 

0.1059 

Temperature  gas  inlet  [K] 

973.15 

Pressure  gas  inlet  [bar] 

1.04 

Fuel  inlet  flow  [mol  s-1] 

2.457  x  10~6 

Fuel  molar  fraction  (equal  to  a 

0.17  CH4,  0.48  H20 

pre-reforming  rate  of  30%) 

0.26  H2.  0.02  CO,  0.04 
C02 

Fuel  steam/carbon  ratio 

2.85 

Fuel/air  flow  time  constant  [s] 

2/- 

Air  inlet  flow  [mol  s-1  ] 

8.282  x  10~5 

Air  molar  fraction 

0.21  02,  0.79  N2 

Average  cell  voltage  [V] 

0.7  (cascaded  cell:  7  V) 

Fuel  utilization  [%] 

85% 

flow  direction  as  a  consequence  of  low  temperatures  and  high 
polarization  at  the  fuel  inlet  of  the  cell.  Even  if  the  temper¬ 
ature  increases  up  to  a  cell  length  of  7  cm,  the  local  current 
density  decreases  from  a  cell  length  of  4  cm.  This  is  because 
the  decreasing  of  the  Nernst  voltage  and  the  increasing  of 
concentration-sensitive  losses  are  larger  than  the  influence  of 
temperature-sensitive  losses.  At  a  cell  length  of  7  cm  the  tem¬ 
perature  reaches  a  maximum  of  1268  K.  The  temperature  drop 
in  the  last  third  of  the  cell  length  occurs  due  to  radiative  and 
convective  heat  transfer  between  air  preheater  tube  and  cathode. 

The  local  cell  voltage  shows  a  slightly  maximum  in  the  mid¬ 
dle  of  the  cell  length,  which  can  be  explained  by  ohmic  losses 
in  the  anode  and  cathode.  This  is  because  the  current  collec¬ 
tors  (Fig.  1)  leads  the  highest  current  density  at  the  beginning 
of  the  anode  (cell  length:  0  cm)  and  at  the  end  of  cathode  (cell 
length:  10  cm).  The  average  ohmic  resistances 

8RS 

^An(Ca)  =  (39) 

of  each  electrode  and  each  cell  slice  are  adjusted  to 
4.23  x  10-4  !T2  to  reach  the  same  average  power  density  as  the 
cascaded  micro-tubular  cell.  This  is  much  lower  than  the  average 
ohmic  resistance  of  the  cascaded  concept.  On  these  conditions 


Fig.  9.  Hydrogen  and  methane  fraction  of  the  cascaded  cell. 


the  common  micro-tubular  cell  reaches  a  voltage  of  0.7  V  at 
an  all  over  current  of  2.3  A  which  is  equal  to  an  average  cur¬ 
rent  density  of  244  mA  cm-2  and  an  average  power  density  of 
170  mW  cm-2. 

In  this  study,  the  cascaded  micro-tubular  cell  should  oper¬ 
ate  on  the  same  effective  area,  average  cell  voltage,  gas  inlet 
flow  rates  and  fuel  utilization  as  the  common  micro-tubular  cell. 
To  fulfil  these  conditions  with  10  cascades,  the  all  over  cell 
voltage  must  be  10  times  higher  and  the  current  must  be  10 
times  lower  than  the  values  of  the  common  cell.  Because  the  all 
over  current  of  each  cascade  must  be  the  same,  the  current  den¬ 
sity  is  distributed  over  the  whole  cell  length.  This  explains  why 
the  maximum  of  hydrogen  fraction  is  displaced  to  the  cell  inlet 
(Fig.  9)  and  why  the  decrease  of  the  hydrogen  fraction  shows  a 
linear  characteristic. 

The  average  current  density  (Fig.  10)  related  to  the  cascade 
cell  area  has  the  same  value  as  the  average  current  density  of 
the  common  cell.  The  cascade  voltage  increases  along  the  flow 
direction  as  a  consequence  of  low  temperatures  and  high  polar¬ 
ization  at  the  fuel  inlet  of  the  cell.  From  a  cell  length  of  4  cm  the 
cascade  voltage  decreases  due  to  the  decreasing  Nernst  voltage 
and  the  increasing  of  polarization  losses. 

The  temperature  profile  of  the  cascaded  micro-tubular  cell 
shows  the  similar  tendency  with  the  common  micro-tubular 
cell.  The  temperature  maximum  of  1266  K  is  located  at  7  cm 
of  cell  length.  On  these  conditions  the  cascaded  micro-tubular 
cell  reaches  a  voltage  of  7  V  at  an  all  over  current  of  0.23  A 


Fig.  10.  Local  voltage  and  current  density  of  the  cascaded  cell. 
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Fig.  1 1 .  ////-characteristics  of  the  common  and  cascaded  cell. 
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Fig.  12.  Load  change  of  the  common  cell. 


which  is  equal  to  an  average  current  density  of  244  mA  cm“- 
per  cascade  and  an  average  power  density  of  170  mW  cm-2  per 
cascade. 

Fig.  1 1  shows  the  average  f///-characteristic  of  the  common 
and  cascaded  micro-tubular  cell.  As  mentioned  before  the  ohmic 
resistance  (anode  and  cathode)  of  the  common  cell  is  adjusted 
much  lower  than  the  ohmic  resistance  of  the  cascaded  cell. 
Hence,  the  average  cell  voltage  (0.7  V)  of  both  cells  is  the  same 
at  an  average  current  density  of  244  mA  cm-2.  At  lower  current 
densities  than  100  mA  cm-2  the  voltage  of  the  common  cell  is 
limited  by  the  Nernst  voltage  at  the  fuel  outlet.  To  guarantee  the 
operation  of  the  common  cell  in  this  performance  range  the  fuel 
utilization  has  to  be  decreased  which  involves  a  lower  efficiency. 

As  shown  in  Fig.  11,  the  cascaded  concept  gives  the  oppor¬ 
tunity  to  operate  at  lower  part  load  conditions  than  the  common 
concept  without  decreasing  the  fuel  utilization.  At  current  den¬ 
sities  higher  than  100mA  cm-2  both  cells  show  similar  perfor¬ 
mance  for  the  ohmic  and  electrochemical  assumptions  of  this 
study. 

3.2.  Transient  behavior 

In  common  energy  supply  applications  electrical  current  and 
voltage  are  controlled  to  react  on  changes  in  load  demand.  In  fuel 
cell  applications  inverter  are  used  to  fulfil  these  requirements. 
To  decouple  the  transient  behavior  of  the  SOFC  from  that  of  the 
inverter,  the  load  change  is  simulated  for  an  increased  current  of 
+10%.  The  transient  electrical  response  and  the  change  in  tem¬ 
perature  profile  are  simulated.  The  inlet  gas  compositions  and 
inlet  temperatures  of  fuel  and  oxidant  are  kept  constant  in  this 
simulation.  In  practical  SOFC  systems  fuel  flow  rates  are  con¬ 
trolled  to  keep  the  steady  state  fuel  utilization  ratio  constant.  This 
implies  that  reacting  species  are  delivered  with  a  time  delay  (Eq. 
(38)).  In  this  study,  the  air  ratio  is  kept  constant.  It  is  assumed 
that  the  time  constant  of  the  fuel  cell  capacities  were  negligible 
in  comparison  to  the  time  constant  of  the  gas  deliveries.  The  per¬ 
fusion  of  the  gas  streams  through  the  cell  occurs  in  a  time  lower 
than  0.1  s.  Therefore,  the  transient  behavior  of  gas  distributions 
along  the  cell  length  was  neglected  as  well. 

As  shown  in  Fig.  12,  the  fuel  utilization  of  the  common  micro¬ 
tubular  cell  increases  up  to  94%  at  a  time  of  0  s.  The  minimum 
cell  voltage  during  this  load  change  is  0.656  V.  The  polarization 


losses  at  the  time  of  0  s  are  much  higher  compared  to  the  losses 
at  the  new  steady  state  condition  because  of  the  temporary  high 
fuel  utilization.  This  is  because  the  cell  voltage  increases  with 
the  decreasing  fuel  utilization.  At  the  new  steady  state  condition 
the  fuel  utilization  is  controlled  to  85%.  The  voltage  reaches  a 
final  value  of  0.687  V.  The  fuel  utilization  could  have  a  limiting 
influence  in  load  change,  if  the  reacting  species  are  delivered 
with  high  time  delays.  This  effect  applies  especially  to  the  com¬ 
mon  micro-tubular  cell,  which  is  from  a  theoretical  point  of  view 
unable  to  reach  cell  voltages  above  the  lowest  Nernst  voltage. 

Fig.  13  shows  the  transient  temperature  profiles  of  the  com¬ 
mon  micro-tubular  cell.  At  a  time  of  300  s  the  temperature 
profile  reaches  90%  of  the  final  profile  (r  =  oo).  The  temper¬ 
ature  maximum  increases  from  1268  K  up  to  1288  K.  Due  to  a 
local  displacement  in  current  density  the  temperature  maximum 
is  displaced  slightly  to  the  cell  inlet. 

As  shown  in  Fig.  12,  the  influence  of  changing  in  cell  tem¬ 
perature  on  polarization  losses  can  be  neglected.  The  outlet 
temperature  of  gas  species  increases  from  1131  to  1142K  at 
the  cathode  and  from  1173  to  1188Kat  the  anode.  In  deviation 
to  this  study  the  air  ratio  would  be  controlled  in  practical  SOFC 
stacks  to  keep  the  maximum  temperature  of  the  cell  constant. 

The  transient  response  of  the  cascaded  micro-tubular  cell 
shows  a  similar  behavior  with  the  common  micro-tubular  cell. 
Even  if  the  temporary  high  fuel  utilization  occurs  during  this 
load  change  as  well,  the  minimum  average  cell  voltage  (0.67  V) 
is  higher  than  the  voltage  of  the  common  micro-tubular  cell 
(0.656  V).  This  can  be  explained  by  the  fact  that  the  change  in 


Fig.  13.  Temperature  profile  of  the  common  cell. 
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Fig.  14.  Load  change  of  the  cascaded  cell. 


Fig.  15.  Temperature  profile  of  the  cascaded  cell. 


polarization  is  mainly  located  in  the  last  half  of  the  cascaded  cell. 
At  the  new  steady  state  condition  the  fuel  utilization  is  controlled 
to  85%.  The  voltage  reaches  a  final  value  of  6.85  V  (Fig.  14). 

The  transient  temperature  profiles  of  the  cascaded  micro¬ 
tubular  cell  are  shown  in  Fig.  15.  At  a  time  of  300  s  the  tem¬ 
perature  profile  reaches  95%  of  the  final  profile  (r  =  oo).  The 
temperature  maximum  increases  from  1266  K  up  to  1282  K. 
Due  to  the  regular  distribution  in  current  density  the  tempera¬ 
ture  profile  shows  no  displaced  at  new  steady  state  condition. 
The  outlet  temperature  of  gas  species  increases  from  1149 
to  1160K  at  the  cathode  and  from  1196  to  1210  K  at  the 
anode. 

In  comparison  to  the  transient  response  of  the  common  micro¬ 
tubular  cell  the  cascaded  concept  shows  a  more  favorable  behav¬ 
ior  with  respect  to  voltage  drop  and  temperature  distribution. 

4.  Conclusion 

A  two-dimensional  model  was  developed  to  simulate  the 
steady  state  and  transient  characteristics  of  tubular  SOFC. 
Steady  state  and  transient  characteristics  of  a  common  and  cas¬ 
caded  micro-tubular  cell  with  a  cell  diameter  of  2  mm  were 
simulated  and  compared. 

In  this  study,  the  cascaded  micro-tubular  cell  has  the  same 
effective  area,  average  cell  voltage,  gas  inlet  flow  rates  and  fuel 
utilization  as  the  common  micro-tubular  cell.  With  10  cascades 
the  all  over  cell  voltage  must  be  10  times  higher  and  the  cur¬ 
rent  must  be  10  times  lower  than  the  values  of  the  common  cell. 


To  reach  a  common  average  power  density  at  an  average  volt¬ 
age  of  0.7  V  the  average  ohmic  resistance  (anode/cathode)  of 
the  common  micro-tubular  cell  were  adjusted  to  4.23  x  10-4  Q. 
This  is  84  times  lower  than  the  average  ohmic  resistance  of 
the  cascaded  concept.  This  resistance  depends  strongly  on  the 
thickness  of  anode  and  cathode,  wire  thickness  and  contact 
resistances.  Even  if  all  over  current  is  reduced  by  the  fac¬ 
tor  of  cascades  and  power  losses  depends  on  all  over  current 
approximately  to  square,  the  current  density  in  ^-direction  could 
reach  very  high  values,  with  a  small  thickness  of  the  anode 
or  cathode.  Compared  to  the  ohmic  resistance  of  the  tubular 
cell  developed  by  Westinghouse  the  adjusted  ohmic  resistance 
of  the  common  micro-tubular  cell  is  this  still  a  low  value. 
However,  future  development  in  cell  concepts  may  reach  the 
adjusted  resistance  of  the  common  micro-tubular  cell  of  this 
study.  In  that  case,  the  common  micro-tubular  cell  could  reach 
the  same  power  density  as  the  cascaded  micro-tubular  cell  of  this 
study. 

The  (///-characteristic  of  the  cascaded  cell  shows  that  higher 
average  voltages  can  be  reached  at  low  average  current  densities 
compared  to  the  common  micro-tubular  concept. 

The  transient  electrical  response  and  the  change  in  temper¬ 
ature  profile  were  simulated.  The  inlet  gas  compositions  and 
inlet  temperatures  of  gas  species  were  kept  constant.  The  fuel 
flow  rates  were  controlled  to  keep  the  steady  state  fuel  uti¬ 
lization  ratio  constant.  The  results  of  the  transient  simulation 
show  the  quick  response  to  a  load  change  in  the  common 
micro-tubular  cell  and  cascaded  concept  as  well.  During  this 
load  change,  the  minimum  average  cell  voltages  reached  val¬ 
ues  of  0.67  V  for  the  cascaded  cell  concept  and  0.656  V  for  the 
common  cell  concept.  It  is  supposed  that  the  cascaded  cell  con¬ 
cept  could  work  on  higher  load  changes  than  the  common  cell 
concept. 

Future  work  should  involve  the  evaluation  of  the  simulation 
results  with  experimental  data.  It  is  assumed  that  the  cascaded 
cell  concept  offers  a  high  potential  in  optimization  with  respect 
to  the  number  of  cascades,  the  anode  and  cathode  thickness 
respectively.  Not  at  least  economic  analysis  including  the  periph¬ 
eral  system  costs  should  be  taken  into  account. 
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